Signature candidate of quantum chaos far from the semiclassical regime 
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We numerically investigated the entanglement product in the simplest coupled kicked top model 
with the spin j — 1. Different from the dynamical pattern of entanglement in the semiclassical 
regime, two similar initial states may have discordant entanglement oscillations. A candidate of 
the quantum signature of this classical chaotic system was proposed. The potential antimonotonic 
relation between the rank correlation coefficient qualifying the concordant of two entanglement 
evolutions and the stationary entanglement was preliminarily revealed. 
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Entanglement in the quantum chaotic systems has 
been studied by several authors [l|-|24| . For the system of 
weak coupled kicked tops, in the semiclassical regime, it 
has been clarified that two initially separable subsystems 
can get entangled in a nearly linear rate depending on 
the intrinsic chaotic properties, and their entanglement 
eventually reach a saturation |8l— fill. |20| . It has also been 
elucidated that the increment of the nonlinear parame- 
ter of weakly coupled kicked tops does not accelerate the 
entanglement production in the strongly chaotic region 
fllj ] . For strong coupled kicked tops in the semiclassi- 
cal regime, it is observed the greater the kick strength, 
the higher asymptotic value of entanglement if the kick 
strength is not very small [10 . |20j. In the deep quan- 
tum regime, based on the combined electronic and nu- 
clear spin of a single atom, the experimental realization 
of the quantum kicked top with J = 3 has been presented 
and it has been found good correspondence between the 
quantum dynamics and classical phase space structures 
[3|. Clear differences in the sensitivity to perturbation 
in chaotic versus regular regimes, and experimental evi- 
dence for dynamical entanglement as a signature of chaos 
have also been observed |3| . In Ref . [25| , it has been shown 
that the discrepancy between quantum and classical evo- 
lutions can be significantly decreased by even a weak loss 
of the coherence. Motivated by previous works, we in- 
vestigated the entanglement quantified by log-negativity 
in coupled kicked tops with the spin j = 1 by using 
^ = p|^i)(^i|®|^2)(^ 2 | + (l-p)|^ 2 >(^ 2 |®|V'i)<V'i|aBthe 
initial state. The difference of the Von- Neumann entropy 
between p p= o and p P =o.5 is assumed to be AS, which de- 
pends on the overlap of \ipi) and {ipz). The purpose of 
choosing this kind of initial states is to keep the energy in- 
dependent on the initial mixedness for the following cou- 
pled symmetric kicked tops. Moreover, this kind of initial 
states guarantee the mixedness-independent dynamical 
trajectories of arbitrary symmetric collective observables 
such as J% +Jz 2 or JziJz 2 m two coupled symmetric tops. 
It is rationally conjectured the characteristics of the en- 
tanglement evolution E 2 (t) corresponding to the above 
initial state with p = 0.5 reflects closer details concerning 



FIG. 1: The log-negativity of two tops is plotted as the func- 
tion of t for three different nonlinear coefficients and for pure 
or mixed spin coherent state, j — 1; e = 0.05; 71 = —3; 
72 = 3; (a) k = 0.25, p = 0.5; (b) k = 0.25, p = 0.0; (c) k = 3, 
p = 0.5; (d) k = 3, p = 0.0; (e) k = 6, p = 0.5; (f) k = 6, 
p = 0.0. 



its classical counterpart than the entanglement evolution 
E\{t) with initial p — 0. For investigating how the ini- 
tially lost information affects the entanglement dynami- 
cal behaviors of two tops with small spin, it is proposed 
to apply the scaled rank correlation coefficient 



lim 
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[Q(E 1 (t)-E 1 (t~l))Q(E 2 (t)-E 2 (t-l))} 

(1) 



to quantify the concordant or discordant between two en- 
tanglement evolution, where @(x > 0) = 1, ©(a; = 0) = 
and <d(x < 0) = — 1. The T^ closer to 1 indicates two 
entanglement evolutions are more concordant. For the 
weak coupled kicked tops with j > 50, previous studies 
have shown the entanglement between two tops does not 
exhibit large amplitude oscillation but monotonic ascent 
to a saturation [8h11|. [201] . In these cases, T^ is always 
near one which can not provide the explicit quantum sig- 
nature of its classical chaotic counterpart. 

The Hamiltonian of the coupled kicked tops can be 
written as M- 
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ft(t) = #i(i) + H 2 (t) + H 12 (t) 



(2) 




FIG. 2: The same parameters as Fig.l; (a) k — 0.25; (b) 
k = 3; (c) k = 6; (Solid line) p = 0; (Dash line) p = 0.5. The 
scaled rank correlation coefficients Yd corresponding to k have 
the order relation T d (k = 0.25) < T d {k = 6) < T d (k = 3). 




FIG. 3: The log-negativity of two tops is plotted as the func- 
tion of t for different nonlinear coefficients and pure initial 
spin coherent state, j — 1; e = 0.05; 71 = —3; 72 = 3; 
p = 0.0; Im(fc) = 0.01; (Solid line) Re(fc) = 3.0; (Dash line) 
Re(fc) = 6.0; (Dot line) Re(fe) = 0.25. 



with H^t) = 1 4 + |i Jl E S(t - n), (3) 



2.7 



#12 (t) = r44E^ 



(4) 





Re[k] 



FIG. 4: The log-negativity of two tops at the kick step 
t — 1000 is plotted as the function of the real part Re(fc) 
of k for two different pure initial spin coherent state, j — 1; 
e — 0.05; p — 0; Im(fc) = 0.01; (hole square) 71 = 72 = 
exp(0.63i) tan(0.89/2), a chaotic point in phase space of the 
single kicked top with k — 3 in Ref. [g|; (Real square) 71 = —3; 
72 = 3. 



where i = 1,2, and Hi(t)'s are the Hamiltonians of the 
individual tops. H 12(f) represents the coupling between 
two tops with a coupling strength of e/j. All these spin 
operators obey standard commutation relations. Choos- 
ing the simultaneous eigenstates \j\,mi)x ® |j2, "7,2)2 of 
the four mutually commuting operators Jf and J Zi (i — 
1,2) as the basis, where Jf\ji,rm)i = ji(ji + l)\ji,mi)i 
and J Zi \ji, triiji = mi\ji,mi)i. Thereafter, it is assumed 
jl = j'2 = j and ki — k2 = k for the coupled kicked 
symmetrical tops. 

In previous studies, usually the pure direct product 
state of two spin coherent state |§-flll [201 ] is chosen as the 
initial state of the total system. Here, we use a mixed 
state of two direct product of spin coherent state as our 
initial state for two tops. The spin coherent state is given 
in \j, m,i)i basis as : 



l7)i = (l + l7lV £ 



v 3— ™« 




where 7 = exp(i(/>) tan(0/2). 
mixed state can be written as 



.(5) 
The initial normalized 



p(0)=p|^ + (0))(^+(0)| + (l-p)|^-(0))(^-(0)|, (6) 

where \ip + (0)) and \ip~(0)) are given by: 

l^ + (0)) = |7i)i|72) 2 

|V"(0)) = l72)i|7i>2, (7) 

The log-negativity of a bipartite density matrix is de- 
fined by [26| 



Af(p)=log 2 \\p T \\ 



(8) 



where p is the partial transpose of p and \\p || denotes 
the trace norm of p , which is the sum of the singu- 
lar values of p r . The coupled kicked tops with j = 1 
is far from the semiclassical regime. Due to the ex- 
change symmetry of tl(t) and its Floquet operator U, 
the Floquet eigenstates have well-defined exchange par- 
ity: j e = (j+l)(2j+l) Floquet eigenstates jej with even 
exchange parity and j — j(2j + 1) Floquet eigenstates 
\oi) with odd exchange parity. The Floquet operator U 
can be expressed as 



U = Y,^'\e l )(e l 



je+3, 

E - 

l=je + l 



t(f)i 



Oi)(Oi 



(9) 



The initial state can be expanded by the Floquet states 
as 



Je 



Je+Jo 



|7'2>i|ti>2 =E a »l e *) + E bl \° l 

i—l i— j e + l 



and 



Je +jo 



l7l)l|72)2 



E a d e *)- E bi \° i 

i—l i=j 6 -\-X 



(10) 



(11) 



where a, and bi are the complex coefficients. The evolv- 
ing state corresponding to the initial states p p= o(0) and 
Pp=o.5(0) can be obtained as 



je +Jo 



n72)i|7i>2 = J> em0, l e *>+ E b ^ m4H \°i)^ (12) 



j=l 



and 



Je+Jo 

+ ^ fti&re^^-^'hXo, 1,(13) 

respectively. Comparing the above two evolving states, it 
is explicit that the scaled rank correlation coefficient T^ 
reflects the influence of coherent superposition between 
two Hilbert subspaces with different exchange parity on 
the dynamical behaviors of entanglement. In Figs. (1-2), 
we can see the log- negativity of the time evolution density 
operator for two different mixedness and for k — 6, k = 3 
and k = 0.25, which are strong chaotic, weak chaotic and 
regular cases of the single kicked top in previous litera- 
tures [27|, respectively. For all of the cases, we pick an ini- 
tial state in Eqs.(6-7) with 71 = —3 and 72 = 3. Slightly 
different to the results in Ref. [3j , the log- negativity shows 
quasi-periodic behavior and collapses and revivals of the 
oscillation for k = 3. While for the case of k = 6, the 
global dynamical pattern of the log-negativity is quasi- 
period but local dynamical pattern is chaotic. For small 
value oik = 0.25, i.e. the regular case, the log-negativity 
only exhibits a chaotic oscillation. For the mixed ini- 
tial state in Eqs.(6-7), the evolving log- negativity with 
p = 0.5 is the lower bound of the one with p = 0. The 
scaled rank correlation coefficient T^ quantifying the con- 
cordant or discordant of the entanglement evolutions cor- 
responding to initial states with p — and p = 0.5 de- 
pends on k even though the corresponding evolving states 
contain the same dynamics for arbitrary symmetric col- 
lective operators such as J^ + J% or the coupling term 

J Zl J Z2 . From Fig. 2(b), it can be found that the corre- 
sponding entanglement dynamical behaviors of two cases 
with p — and p — 0.5 are synchronous and ordered 
when k = 3, not like the disorder and asychronous corre- 
spondence in Fig. 2(a) and Fig. 2(c). In Ref.[3j, it has been 
experimentally demonstrated that less average entangle- 
ment generation for initial states localized in regular re- 
gions compared to those in the chaotic sea. Similarly, 
taking into account of the external pulsed monitoring the 
collective observable J% + J^ 2 of two spins which makes 
the states of two spins collapse toward the subspace with 
i{J% ( )i = j 2 (i = 1,2, and i(Y)i denotes the expectation 
value of the operator Y), we investigated the entangle- 
ment production in the coupled nonhermitian kicked tops 
by adding a positive imaginary part of k in the Hamilto- 
nian in Eq.(3). The results are shown in Figs. (3-4). Both 



the log-negativity of the tops with k = 0.25 and k = 6 
exhibit the damped chaotic oscillation and converge to 
be stationary. The stationary log-negativity firstly de- 
creases with the increase of Re(fc) from to 3, and reaches 
a fiat bottom, then increases again and achieves a local 
maximal value near 6.4. For further increase of Re(fc), 
it is conjectured similar cycles repeat. From Fig. 4, it 
can also be observed the stationary log-negativity is not 
very sensitive to the initial state. Comparing the Fig. 2 
and Fig. 4, it is seen the larger (smaller) stationary log- 
negativity, the more chaotic (regular) dynamical behav- 
iors of the log-negativity, conforming to the conclusion 
in Ref. [J-S [2I| ■ Thus it is proposed to use this station- 
ary log-negativity or T^ defined in Eq.(l) as a candi- 
date of the signature of quantum chaos in deep quantum 
regime. Obviously, the strength of the imaginary part 
Im(fc) will affect the stationary log-negativity. Usually, 
for a fixed value of Re(fc) and other specific parameters, 
the stationary log-negativity versus Im(fc) exhibits a few 
period of oscillation firstly, then slightly decreases with 
a very small slope, and finally freezes. If the stationary 
log-negativity was adopted as the signature to investigate 
the border between the regular and chaotic regions in the 
stroboscopic quantum phase space or parameter space, it 
is suggested to choose a small enough imaginary part of 
k, i.e. Im(fc) <C Rc(fc) for avoiding to significantly disturb 
the system of interest. 

In Fig. 2, it can be observed the log-negativity of two 
tops achieves a local maximal value near t = 70 for dif- 
ferent values of the parameters k and p. We depict the 
dependence of the local maximal log-negativity on the pa- 
rameter p for different values of k in Fig. 5. They decline 
with p at the slope of about 0.9 for the case with 71 = —3 
and 72 = 3. The slopes only slightly depend on the pa- 
rameter k, but heavily correlate with the initial points 
labeled by 71 and 72. In the cases of 71 = —72, Figs. (6- 
7) show the slopes decrease with the increase of absolute 
value of 71 due to the mixedness of the initial state de- 
creases with I71I for p £ (0, 0.5]. For those kind of initial 
states in Eqs.(6-7), there may be quantum correlation be- 
tween two tops though they are separable |28l . |29| . Their 
quantum discord is dependent of the overlap (71I72) and 
p. The intermediate value of the overlap and p could 
make two tops achieve the larger quantum discord. It 
is also interesting to investigate how the quantum dis- 
cord and classical correlation evolve in the kicked cou- 
pled tops, which maybe embody good quantum-classical 
correspondence [30l - [33l |. 

At last but not least, we investigate the fidelity decay 
in the coupled kicked tops. Usually, the fidelity for two 
mixed states p\ and P2 is defined by [34| 



F(pi,p2) = [Tr(J y/pTp 2 y/pi)] 



(14) 



Here, we adopt this fidelity to study the sensitivity of 
evolving state in the coupled kicked tops to the small 
perturbation, i.e. pi(t) and P2(t) represent two evolv- 
ing density matrix of the kicked coupled tops with two 
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FIG. 5: The dependence of log-negativity at the kick step 
t — 70 on the parameter p is depicted for different values of 
k. The initial states in Eqs.(6-7) with values of 71 = —3 and 
72 = 3 is chosen as an illustration, j — 1; e = 0.05. 
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FIG. 6: The dependence of log-negativity at the kick step 
t — 70 on the parameter p is depicted for the initial states in 
Eqs.(6-7) with different values of 71 and 72. j — 1; k — 3; 
e = 0.05. 




FIG. 8: The fidelity of the system is plotted as the function 
of t with j = f , k = 0.25, k' = 0.26, 71 = -3, 72 = 3. (a) 
e/j = 0.02, (Solid line) p = 0; (Dash line) p = 0.2; (Dot 
line) p = 0.5. (b) e/j = 0.05, (Solid line) p = 0; (Dash line) 
p = 0.2; (Dot line) p = 0.5. 
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FIG. 7: The dependence of log-negativity at the kick step 
t — 70 on the parameter p is depicted for the initial states in 
Eqs.(6-7) with different values of 71 and 72. j — 1; k — 3; 
e = 0.05. 



FIG. 9: The fidelity of the system is plotted as the function 
of t with j = 1, k = 0.25, k' = 0.26, 71 = -3, 72 = 3. (a) 
e/j = 0, (Solid line) p = 0; (Dash line) p = 0.5. (b) e/j = 0.5, 
(Solid line) p = 0; (Dash line) p = 0.5. 



slight different nonlinear coefficients k and k' = k + Sk 
{5k -C 1), respectively. Initially pi(0) = p2(0) belong to 
those class of states in Eqs.(6-7). Generally, the fidelity 
decay should be taken average of the overall initial states 
in the phase space, and the average fidelity decay rate 
can be built relations with the signatures of its classical 
chaotic counterpart such as the mean Lyapunov exponent 
35 37, 39]. The average fidelity decay typically exhibits 
three consecutive stages, i.e. the short-time parabolic 
decay, the followed intermediate-time asymptotic decay, 
and long-time saturation |36l-|38|. Nevertheless, it is still 
worthy to analyze the corresponding fidelity decay for 
some specific evolving states with manifested chaotic dy- 
namical behaviors of entanglement. The key parame- 
ters in Fig. (2a) are adopted to reveal how the fidelity 
decay in such a case with chaotic dynamical behavior 
of bipartite entanglement. The results are displayed in 
Figs. (8-9). In the first stage, the fidelity more rapidly 
decays than exponential decay, and the decay rate de- 
creases with the increase of mixedness and the coupling 
coefficient e. Then, the fidelity revivals. When e — 0, it 
can be observed perfect recurrence of the fidelity, which 
implies the uncoupled kicked top with k = 0.25 is regu- 
lar in the long-time scale. In the strong coupling case 
with e = 0.5, initially the two trajectories of fidelity 
decay seem coincident, but after some time, the diver- 
gence is obvious and enlarged. In the kicked coupled 
spin-1, the coupling coefficient e plays the more crucial 
role in the chaotic dynamical behaviors of the entangle- 



ment than the nonlinear coefficient k. Via replacing the 
operators of the spin j = 1 by Ji = J2 



ailjl and 



Jy — E n =i a rh*/2i where a z „ and Uy„ are the Pauli op- 
erators of the qubit with label {i = l,2;n = 1,2}, the 
kicked coupled spin 1 is equivalent to a kicked four-qubit 
2D Ising square. Four-qubit 2D Ising square may be one 
of the basic cells for constructing the universal quantum 
gate or measurement-based quantum computation pfj| . 
Therefore, to fully investigate the border between the 
regular and chaotic regions both in the quantum phase 
space or parameter space is very important and desirable 
in the future work. 

In summary, this Letter deals with the dynamical pat- 
tern of entanglement in the simplest kicked coupled sym- 
metric tops with j = 1. The main result is to prelim- 
inarily reveal the potential antimonotonic relation be- 
tween the scaled rank correlation coefficient T^ defined in 
Eq.(l) and the stationary entanglement induced by suit- 
able nonhermitian term in the coupled kicked tops. How- 
ever, further calculations including many enough initial 
states and key parameters in the system of the coupled 
kicked tops are necessary to finally qualify these signature 
candidate of quantum chaos in low dimensional Hilbert 
space. Meanwhile, there is still lack of a reasonable phys- 
ical interpretation why the rank correlation coefficient 
has such an antimonotonic relation with the stationary 
entanglement in the system of amplified kicked coupled 
spin-1. 
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